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Abstract 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning 
tunneling microscope (CryoSTM) is proposed. TIA with the tip-sample component 
in CryoSTM is called as CryoSTM-TIA. With transimpedance gain of 1 GQ, the 
bandwidth of the CryoSTM-TIA is larger than 200 kHz. The distinctive feature of 
the proposed CryoSTM-TIA is that its pre-amplifier is made of a single cryogenic high 
electron mobility transistor (HEMT), so the apparatus equivalent input noise current 
power spectral density at 100 kHz is lower than 6 (fA)?/Hz. In addition, “bias-cooling 
method” can be used to in-situ control the density of the frozen DX- centers in the 
HEMT doping area, changing its structure to reduce the device noises. With this 
apparatus, fast scanning tunneling spectra measurements with high-energy-resolution 
are capable to be performed. And, it is capable to measure scanning tunneling shot 
noise spectra (STSNS) at the atomic scale for various quantum systems, even if the 
shot noise is very low. It provides a powerful tool to investigate novel quantum states 
by measuring STSNS, such as detecting the existence of Majorana bound states in the 
topological quantum systems. 


1 Introduction 


High performance transimpedance amplifier (TIA) for cryogenic scanning tunneling mi- 
croscope (CryoSTM) is a key element [I2]. In Ref.[3], it was proposed a design of TIA for 
CryoSTM. TIA with the tip-sample component in CryoSTM is called as CryoSTM-TIA. 
In that design, the CryoSTM-TIA has a transimpedance gain of 1 GO, a bandwidth of 
more than 300 kHz, and its equivalent input noise current power spectral density (PSD) 
that characterizes the circuit inherent noise is 21 (fA)?/Hz at 100 kHz. The CryoSTM- 
TIA is capable of measuring the scanning tunneling shot noise spectra (STSNS) of various 
quantum systems at atomic resolution, even if the measured noise current PSD is as small 
as a few (fA)?/Hz. In Ref.[3], it was illustrated how the CryoSTM-TIA is used to mea- 
sure STSNS of quantum systems to investigate novel quantum states, such as detecting 
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Majorana bound states (MBSs) in iron-based superconductors [4]. In the CryoSTM-TIA, 
cascade the pre-amplifier (Pre-Amp) and the post-amplifier (Post-Amp) to form an oper- 
ational amplifier (OPA) denoted as Macro-OPA. Pre-Amp in Macro-OPA is a differential 
amplifier made of a pair of high electron mobility transistors (HEMTs), where the HEMT 
is a kind of cryogenic GaAs MESFET [51 [6]. However, it is quite difficult to select cryogenic 
GaAs MESFETs with identical performances for pairing [2]. 

In this work, a design for the CryoSTM-TIA with a single HEMT in Pre-Amp is 
proposed. In the CryoSTM-TIA, cascade Pre-Amp and Post-Amp to form an inverting- 
amplifier (Inv-Amp), which is not OPA. The CryoSTM-TIA still has transimpedance gain 
of 1GQ and bandwidth more than 200 kHz. However, its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz, i.e. 1/4 of that in Ref.[B], since a noisy HEMT is 
reduced. Therefore, the lower tunneling shot noise of quantum systems can be measured 
at the atomic scale with this apparatus. Furthermore, only a single HEMT is used in 
Pre-Amp, avoiding the difficulty of matching two identical HEMTs. And, “bias-cooling 
method” can be used to in-situ reduce the inherent noise of the HEMT. With these 
advantages, this apparatus will be a powerful tool to investigate novel quantum states in 
various quantum systems, in the wider range of applications than that in Ref.|3], such as 
high-T, superconductors [7 [8], topological superconductors [4], MBSs [4 {9} [12], 
Andreev reflection [13} [4], and Kondo effect [15], etc. 


2 Circuit of the proposed CryoSTM-TIA 


The circuit of the proposed CryoSTM-TIA is shown in Fig[i] It consists of several compo- 
nents: the single HEMT amplifier part of Pre-Amp shown in the dashed box (a1) of Figf1| 
the power supply of Pre-Amp in the dashed box (a2), Post-Amp shown in the dashed 
box (b), the compensated feedback network shown in the dashed box (c), and the signal 
source circuit shown in the dashed box (d). The two stage amplifier made of Pre-Amp 
and Post-Amp is called as Inv-Amp. Inv-Amp is connected with the feedback network 
to form TIA. The input G of TIA is connected to the signal source circuit, to form the 
CryoSTM-TIA. The components placed in the cryogenic zone are shown in the dotted 
box. The parameters of all components of the CryoSTM-TIA circuit are listed in Table [I] 


2.1 Design of Pre-Amp 


As shown in Fig[i] only a single cryogenic HEMT is used in Pre-Amp. The HEMT in 
this work is CNRS-HEMT (denoted as HEMT) [5] [6] developed by CNRS/LPN in France 
with the excellent cryogenic and noise characteristics. It is capable of operating under 0.5 
K with only 0.1 mW on the ideal operating point of “Vas = 100 mV and Ias = 1 mA”. Its 
parameters are listed in Table j! where e?, is its equivalent input noise voltage PSD and 
is is its equivalent input noise current PSD [6]. 

The single HEMT amplifier part of Pre-Amp is shown in dashed box (al) in Fig/1] 
The source of HEMT is grounded through a resistor Rs and a small variable resistor Rs1. 
Rs = R; + Rg, and a capacitor Cs of 11 uF is in parallel with Rg. Drain O; of HEMT 
is connected to the load resistor R,. The resistor Rp is grounded through the variable 
resistor Rp. Ry = Ry + Ry. The other end of Ry, denoted as O2, is connected to the 
resistor R2, and L as the other end of Rə is connected to P (the other end of Rı) through 
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Figure 1: Circuit of the proposed CryoSTM-TIA. Single HEMT amplifier part of Pre- 
Amp is shown in the dashed box (al), the power supply of Pre-Amp in the dashed box 
(a2), Post-Amp in the dashed box (b), the compensated feedback network in the dashed 
box (c), and the signal source circuit in the dashed box (d). The components placed in 


the cryogenic zone are shown in the dotted box. The parameters of all components of 
CryoSTM-TIA circuit are listed in Table [I] 


a small variable resistor Rp. Ri = R2 = Ry. P is connected to the power supply of Pre- 
Amp. P and L are grounded through the large capacitors C1 and C2 respectively. In the 
power supply, Rp = Ry + Ry, and Ry is a variable resistor. C'a is the input capacitance 
of Pre-Amp. HEMT, R,, Rp, Ri, and Rə are placed in the cryogenic zone, shown in the 
dotted box. C{ is the capacitance of the cable that connects Gate G of HEMT and the tip 
in the CryoSTM. HEMT is as close to the tip as possible to reduce Cy, which can be less 
than 0.5 pF [16]. The operating point adjustment for Pre-Amp is shown in Sect [4] The 
Pre-Amp parameters are shown in Table [I] 

For Pre-Amp, as AC voltage V is applied to input G, the AC voltage difference between 
Oı and Og is denoted as View: Ayp = Vop j V is the voltage gain of Pre-Amp. Since the 
trannsit frequency for the CNRS-HEMT is gm/(27Cg,) © 1 GHz, the bandwidth of Pre- 
Amp is more than 30 MHz. In max{gm/(27Cs), 1/(2rRsCs)} << f < 3 MHz (i.e. 30 Hz 
<< f < 3 MHz with the parameters listed in Table[lp, 


Ayp x —9mfRa, (2.1) 


where Rg = RL/(1 + gaR). ie, Ayp ~ —20 in 30 Hz<< f < 3 MHz. The input 
capacitance of Pre-Amp is 
Ca = Ces + (1 — Ave) Cga- (2.2) 


Table 1: Parameters of all components of CryoSTM-TIA circuit 


HEMT (CNRS-HEMT) 


Gate-—source resistance Ra 
Transconductance gm 
Channel conductance gq 
Gate-source capacitance Cs 
Gate-drain capacitance Cga 
Drain-source voltage Vas 
Drain-source current Ias 


>10 TQ 
40 mS 
1 mS 
5 pF 
1 pF 

100 mV 
1mA 


z 
Equivalent input noise voltge PSD ex n ee aes 
2 
Equivalent input noise current PSD iZ o < a i 
Pre-Amp 
Rs + Rsı 50452 Cs 11 uF 
Ry + Roi 15045 9 Rı, R2 1 kQ 
Rp 5+5 Q 
Ci 22 uF C2 11 pF 
PVR LM4050-10 Rp 2 kQ 
R; + Rs 4400 + 10 Q Vep +15 V 
THS4021 as Rear-OPA 
Gao 94 dB fo 16 kHz 
Ca 1.5 pF Ra 1 MQ 
VRp +15 V Van -15 V 


Equivalent input noise voltge PSD e2 2.25 (nV)?/ 


Hz in f > 10 kHz 


Equivalent input noise current PSD i2 4 (pA)?/Hz in f > 10 kHz 


Post-Amp 

Re 200 kQ Cii, Ci 150 pF 

Rea, Rez 100 Q Cri, Cro 50 pF 
Feedback network 

Rp 1GQ CF ~3 pF 

Ry 300 kQ Ck ~ 0.2 pF 

Ce 10 nF 

Signal source circuit 
Ry >1 MQ Cr ~ 0.5 pF 


Note: + indicates the variable resistance range. Without specification, 


the default value after + is 0. 


i.e., Ca © 26 pF in 30 Hz<< f < 3 MHz. The input resistance of Pre-Amp Ra is the 

gate-source resistance of HEMT, and Ra > 10 TQ, so it can be considered as infinity. 
The power supply of Pre-Amp is shown in the dashed box (a2) in Fig/1] A precision 

voltage reference (PVR) LM4050-10 is used in the power supply to provide a constant 


voltage of 10 V with the typical temperature coefficient of 4 
PVR are eliminated by C1, C2, and Rr. 


-40 ppm/°C. The noises from 


For a differential amplifier, there is a feedback between the two brenches of the transis- 


tors, so the hgih common mode rejection ratio (CMRR) can 


be realized for the differential 


amplifier. For the differential amplifier as Pre-Amp in Ref.|3], if a resistor is instead of 
HEMT (H2) with the same DC resistance, Inv-Amp may become a non-inverting amplifier 
as f — 0. For the design of Pre-Amp in this work, there is no feedback between the HEMT 
brench and Ry brench, therefore Inv-Amp in this work is always an inverting amplifier. 


2.2 Design of Post-Amp and composition of Inv-Amp 


In Figfi] the Post-Amp circuit in the dashed box (b) is the same in structure as that in 
Ref. [B]. There is a commercial operational amplifier (OPA) in the circuit, called as Rear- 
OPA, which is the OPA with high gain-bandwidth-product, such as THS4021, OPA657, 
and LMH6624, etc. The Rear-OPA in this work is THS4021 [18]. The Post-Amp parame- 
ters are listed in Table[i] Ra and C, are the input resistance and capacitance of Rear-OPA 
respectively. Re of 200 kQ is the feedback resistor, and it is placed in the cryogenic zone. 
Post-Amp also contains two cables that connect Q, and Q% in room temperture zone to 
outputs O; and O2 of Pre-Amp in cryogenic zone respectively. The capacitance of the two 
cables is Ciy and Cig respectively, and they may vary from 50 pF to 150 pF. In this work, 
Cii = Cio = C; = 150 pF is assumed. 

A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance of 
100 pF and a distributed inductance of 250 nH, and its distributed LC structure results in 
several resonant gain peaks in [50 MHz, 1 GHz] [19]. Cascade Pre-Amp and Post-Amp to 
form Inv-Amp. In Inv-Amp, Post-Amp and Pre-Amp are connected by the cables (Cable 
O,Q, and Cable 02Q,), and the resonant gain peaks may cause the self-oscillations, as 
some unexpected electromagnetic couplings may exist among the outputs and the inputs 
of the amplifiers. Two low-pass filters (Re = Reg = Re = 100 Q and Cy = Cr2 = Cr = 50 
pF) with the upper cut-off frequency of 30 MHz are added at the Rear-OPA inputs, these 
resonant gain peaks can be greatly reduced, thus avoiding the self-oscillations. 

For the AC signal, the voltage gain of Inv-Amp is aa = aa(f). aa(f) can be expressed 
as 

aa = AypAyr. (2.3) 


By the nodal analysis method, aa can be obtained. And then, with Ap expressed by 
Eq. (2.1), Ayr can be obtained by Eq. (2.3). In max{gm/(27Cg),1/(27RsCs)} << f <3 


MHz (i.e. 30 Hz << f < 3 MHz with the parameters listed in Table[1), 
Rg 1 + j2r f RyiCir , 1 
Fan, TEAG. TI p aE 


aa Ry, 


AvR x 


(2.4) 


where Ry, = Ry RL/(Ry +R), Cir = Ci+C;,, and aa is the voltage gain of the Rear-OPA. 
In this work, Rs >> Ry, and Ra >> Ry. And, ag can be approximately expressed as 
da = aa0o/(1 + j f/f) in (0, 40 MHz] for THS4021. 

With the parameters in Table}1| the Inv-Amp performances are simulated by TINA- 
TI P0]. Fig[] show the curves of |aa(f)|ap and Z (aa(f)) simulated by TINA-TI. By the 
nodal analysis method, aa(f) is calculated by the equations established with all compo- 
nents listed in Table |1| [21], which is also shown in Fig || The curves of |aa(f)|ap and 
Z (aa(f)) calculated by the nodal analysis method are identical with the simulated ones 


respectively. |aa(f)|ap and Z (aa(f)) are also calculated by Eqs.(2.1), (2.3), and (2.4), 


and the calculated curves in Fig [2] are almost identical with the simulated ones respectively 
in [3 kHz, 10 MHz], which verifies the correctness of Eqs. (2.1),(2.3), and (2.4). 


2.3 Frequency compensation of feedback loop 


In order to increase the bandwidth of the CryoSTM-TIA, for the high feedback resistor 
Rp with parasitic capacitance Cp, frequency compensation must be used in the feedback 
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Figure 2: The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|ap, 
and the dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI 
simulation results, the green curves are the calculated results with the nodal analysis 
method, and the black curves are the calculated results obtained by Eqs. (2.1), (2.3), and 


EJ. 


loop [I| 22]. In Figf1| the compensated feedback network is shown in the dashed box 
(c). Taking Ce equalto kCp, where k is above 10°, adjust Ry equalto Rp/k, realizing 
RCs = RpCp. The output voltage of TIA as V, generates the current Ip flowing to input 
G of TIA, so 

Vo Ret Rr Rp 

İp 1+ j2mfRCy 1+ j2rf ROCk’ 


where Cy is the parasitic capacitance of Ry 22]. Zp(f) can be considered as the 
impedance of the feedback network. In Ref.{22], it has been achieved in experiments to 
broaden the bandwidth of the feedback network with the very high feedback resistor Rp 
of 10 GQ to MHz. In (0, 1 MHz], with the parameters listed in Table |1| |Zp(f)| ~ 
Rp/|1 + j2r f RkCr| > Rp/1.008 and |Zp(f)| < Rp, so it can be considered that Zp(f) is 
equal to Rp. 


Zp(f) = 


L pi iiriil L toti tiiid L pouiiit 
2 3.4567 2 2 4567 2 2. 4:35-67 
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Figure 3: The experimental results for the frequency compensation. At 4.2 K, Rp is 1.16 
GQ and Cf is estimated as 3.36 pF. As R, = 380 kQ and Ce = 10 nF, the black curve 
shows that the bandwidth of |Zp/Rplap is broadened to 905 kHz, and the red curve shows 
Z(Zp) as the phase of Zp is —31.5° at 905 kHz. 


In order to verify the above frequency compensation method, the following experiment 


is done. A resistor with 500 MQ(+25 ppm/°C) at room temperature is chosen as Rp, which 
is 1.16 GO at 4.2 K. Rx is a 390 kQ resistor in series with a0 ~ 20 kQ potentiometer, and Ce 
is a COG ceramic capacitor of 10 nF. In order to easily adjust the frequency compensation, 
a 3 pF capacitor is paralleled with the Rp, so that the Cp is about 3.36 pF. Fig.3 shows 
the experimental results for the frequency compensation, in which |Zp| = Rp/V2 at 905 
kHz, and (Zp) as the phase of Zp is —31.5° at 905 kHz. 


2.4 Circuit stability of the proposed CryoSTM-TIA 


In Figfi] the signal source circuit is shown in the dashed box (d), and its parameters are 
shown in Table The differential resistance of the tip-sample tunnel junction (TJ) is 
Ry, which is limited to no less than 107° Rp in this design. The capacitance of TJ is Cy, 
which is in parallel with Ry. And, Cy is estimated as several fF [3]. C = Ca +C,+Cy in 
this work. Cy is at least two orders of magnitude less than Ca + C1, so it can be ignored 
in C and C ~ Ca + Cy. The DC bias & modulated signal voltage source is denoted as 
BMS, which provides the DC bias V; and AC sinusoidal modulated signal voltage V; for 
the CryoSTM-TIA. In the following simulation, Ci is always taken as 0.5 pF [16]. TIA 
connects the signal source circuit to form the CryoSTM-TIA. According to the parameters 
in Table}1} the performances of the CryoSTM-TIA can be simulated by TINA-TI. 
The loop gain Ty, of the proposed CryoSTM-TIA is 


Ti(f) = aa (PEE) = aa (P/A, 


in which 8(f) is the feedback factor, and its reciprocal is, 
1/B(f) ~ 1+ Zpll/Ry +1/Ra + j2rf(Ca + Cy). (2.5) 


By Eq. (2.5), |1/8(f)| and < (1/G(f)) can be calculated. Fig/4{a) shows the calculated 
results of |1/6(f)|ag and < (1/G(f)) with Ry = +00, and Fig/4{b) shows those with 
Ry = 1 MQ. |aalap and Z (aq) are also shown in Fig/4] Both figures show |T.|an = 
laalap — |1/8(f)lap < —10 dB in f > 465 kHz and Z (TŁ) = Z(-aa) — Z(1/8) = 
Z (aa) — 180° — Z (1/8) > —130° in f < 465 kHz. Therefore, the CryoSTM-TIA is stable 
with gain margin more than 10 dB and phase margin more than 50°. 


2.5 Voltage gain and transimpedance gain of the proposed CryoSTM- 
TIA 


With the compensated feedback network mentioned in Sect {2.3} it can be considered that 
Zp is equal to Rp in (0, 1 MHz]. Considering the TJ capacitance Cy, the TJ impedance 
should be Zj = Ry/(1+ j2afRyCy). And, C ~ Ca + Cj. As the AC input voltage V; is 
applied by BMS, the output voltage of the CryoSTM-TIA is Vo, and the voltage gain of 
the CryoSTM-TIA is Ay = Vo/V;. In (0, 1 MHz], by the nodal analysis method, Ay is 


Rp 1 
A% . 3 (2.6) 
y 1 R R ; Rp(Ca+C) 
45 1 aA AF TARA jenf 5 R i 


Setting V; = 0 and applying a sinusoidal current source I, in parallel with TJ, the 
output voltage V, is generated at the output of the CryoSTM-TIA. A; = Vo/J; is called 
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Figure 4: TINA-TI simulation results for the voltage gain Inv-Amp a4, and 1/8(f) (a) with 
Ry = +00, (b) with Ry = 1 MQ. Both figures show |Tilap = |aalap — |1/8(f)lap < —10 
dB in f > 465 kHz and Z (TŁ) = Z(—aa) — < (1/8) = Z (aa) — 180° — Z (1/8) > —130° 
in f < 465 kHz. Hence, the CryoSTM-TIA is stable enough. 


as the transimpedance gain of the CryoSTM-TIA. In (0, 1 MHz], A; is 


Rp 
1 R R . R = atc) 
1 aA aR ARA j2nf 5 7 z 


Disconnecting TIA with the signal source circuit, and sete a sinusoidal current 
source [yp into the input of TIA, the output voltage Vor is generated at the output of 
TIA. Ait = Vor/ Typ is called as the transimpedance gain of TIA. In (0, 1 MHz], A;r is 


R 
Ajr xX E 


JES’ 


Considering Ra >> Rp and |aa| >> 1, Air is 
Rp 


Ait xX 
© 1 jd f RECA’ 
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Figure 5: TINA-TI simulation results for Air/Rp, Air is the transimpedance gain of TIA 
in the proposed CryoSTM-TIA. |Air(fnr)/Rrlag = —3 dB and Z (Ait(fnr)/Rr) = 116° 
at far about 201 kHz. 


In Fig} |Air/Relap and Z (Air/Rp) calculated by Eq. (2.8) with the parameters in 
Table |1} is consistent with their TINA-TI simulation results in [10 mHz, 1 MHz], which 
verifies the correctness of Eq. (2.8) for A;r. The upper cut-off frequency of Air(f)/Rp, 
i.e. its —3 dB frequency fpr is about 201 kHz. Comparing Eq. (2.7) with Eq. (2.8), as 
Ca +Cr+Cy ~ Ca + Cl © Ca and Rj > 107° Rp, the upper cut-off frequency of the 
CryoSTM-TIA is approximately equal to fpr. 


2.6 ‘Transient response of the proposed CryoSTM-TIA 


The simulation results for the transient response of the proposed CryoSTM-TIA are shown 
in Supplemental file 3 [24]. For the CryoSTM-TIA, the time taken from adding the input 
step signal voltage to the output response stably within a certain error is called transient 
response time tr. For the simulations, a resistor with a constant resistance Rg is instead of 
the tip-sample junction. As Ro > 30 MQ, for 0.1% error, tp < 5 us. As 5 MQ < Ro < 30 
MQ, for 1% error, tr < 5 us. As 1 MQ < Ro < 5 MQ, for 3% error, tr < 5 ps. 


3 Inherent noise of the proposed CryoSTM-TIA 


For the circuit of the proposed CryoSTM-TIA shown in Figfi] the differential equivalent 
circuit with all noise sources is used to calculate its equivalent input noise. The details for 
the noise calculations are shown in Supplemental file 4 [25]. 


3.1 Equivalent input voltage noise and equivalent input current noise of 
Inv-Amp 


The equivalent input noise voltage and equivalent input noise current of HEMT are de- 
noted as ey and iy respectively. The resistors Ry, R2, and Rş are in the cryogenic zone 
of 4.2 K, and their noises in f > 1 kHz are thermal noise, which can be neglected [3]. 
Ry is in the cryogenic zone of 4.2 K, and the thermal noise voltage of the resistor Ry 
is erg. The thermal noise voltage of the resistor Re, is e1, and that of Rez is eg. The 
equivalent input noise voltage and equivalent input noise current of the Rear-OPA are 
denoted as eg and ia respectively. These noise sources are independent. The equivalent 


input noise voltage and equivalent input noise current of Inv-Amp are denoted as ea and 
ia respectively. e4 is the equivalent input noise voltage PSD of Inv-Amp, 7% is its equiva- 
lent input noise current PSD, ei, is its equivalent input noise voltage-current PSD, and 
iae% is its equivalent input noise current-voltage PSD. By the nodal analysis method and 
Wiener-Sinchin theorem, ignoring the minor terms, it is obtained that 


a ea ae 2 
B= + Min Stn eee + eS (1 hey (3.1) 
2 T AZ | j Ra 27 P 


E C2 /R 
RoR Puanga) 


A2 R2 
Rut 2 72 
2 a 
+(2rf) (Con + Coa + Ca) g (3.2) 
C R 
Al =e) = jarf e Ez i a++) 
; RyL RyL\ Ë 
= 72 s —— C e 3.3 
jarf (Ca + Cpa + Ra a) 14 2) 4 (3.3) 


According to the parameters listed in Table |1| for Inv-Amp in this work, e2 = = 0.28 
(nV)?/Hz and iå = 0.14 (fA)?/Hz at f = 10 kHz, and “as = 0.09 (nV)?/Hz and a =4,9 
(fA)?/Hz at f = 100 kHz. For the Macro-OPA in Ref. [3], e2 = 0.5 (nV)?/Hz and 74 = 0.8 
(fA)?/Hz at f = 10 kHz, and BA = 0.14 (nV)? /Hz and 2 = = 20 (fA)? /Hz at f = 100 kHz. 
The resistor Ry in Inv-Amp in this work replaces the transistor H2 in the Macro-OPA in 
Ref. [B]. The noise generated by Ry is 2 orders of magnitude smaller than that generated 
by H2 in the Macro-OPA in Ref. [3], i.e. en << A? nee. Therefore, the equivalent input 
noises of Enh in ae ES are much lower than those of the Macro-OPA in Ref. [3]. 

In Egs.(3.1] , and B3), ZyRA / RÌ is one order of magnitude smaller than e. 
In ray a i + m TA is 2 orders of magnitude smaller than e?,. In Eq. (3.2), 
(Cys + Coa + Ca Rut /Ra)* 22/92, is one order of magnitude smaller than C%e2/A2p. In 
Eq. (3.3), (Cos + Cga + Ca Rur /Ra) (1+ Rut/Ra) i2/92, is one order of magnitude smaller 
than Cae?/A2p. Further ignoring the minor terms in Eqs. (3-1), (3.2), and (3.3), 


eh © eh + (e2 + ef + 63)/Ap, (3.4) 
i ws ih + (Qn f) ORE + e? + €2)/A2p, (3.5) 
eaty © (tae4)* = —j2mfCa(e2 + e? + €3)/Arp. (3.6) 
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3.2 Equivalent input voltage noise and equivalent input current noise of 
the proposed TIA 


The equivalent input noise voltage PSD of TIA is denoted as e2, its equivalent input 


noise current PSD is i2, its equivalent input noise voltage-current PSD is epi, and its 
equivalent input noise current-voltage PSD is ipe. They are 


ep =e ~ ef + (e2 +e} + €5)/Abp, (3.7) 


i = 4kpT/Rp + 2B / R2 + 


2 T p2 e 
+i + (27 f)? CK(e2 + a + A (3.8) 


erin, = (ieh) = e} /Rp + eaii 
a+ete 
Ap Rp 
—j2m fC (e? + e? + €2)/Atp. (3.9) 


~ ef /Re + 


3.3 Equivalent input current noise of the proposed CryoSTM-TIA 
The equivalent input noise current PSD of the proposed CryoSTM-TIA is 


2, t=? 4 [1/R5 + (27 f) 202] e2 e, 
+ (1/Ry + j2afCty)eriz + (1/Ry — j2afCs)irer, (3.10) 


where Ciy = C1 + Cy. Putting Eqs.(3.7| , and (3.9) into Eq. (3.10), 


E a a 


+2rf)? |c? Zea + C7(e2 + e? + 6B) /Azp| , (3.11) 
Le. = o 
= iå + 4kBT/Rr + (1/Ry +1/Rp) e + iy, (3.12) 
where a = (2r f}? [oke e + 2C', Ciy (e? + e + e2) /A2p 


For the proposed va ner Rp = 1 GQ, Ca = 26 6 pF, CL = 0.5 Be and Cy = 10 
fF. Rp and TJ are in the cryogenic zone at 4.2 K. As Ry = 1 MQ, i?, and its four 
components are listed in Table [2] The noise components of the CryoSTM-TIA in Ref. [3] 
are also shown in Table The CryoSTM-TIA proposed in this work has an equivalent 
input noise current PSD of 0.62 (fA)?/Hz at 10 kHz and 5.3 (fA)?/Hz at 100 kHz, which 
are much lower than those in Ref.[3}. 

For the Epi-wafer of the CNRS-HEMT, as Si is doped in Al,Ga,_,As, each dopant 


gives rise to two types of electronic states: a shallow and delocalized donor level associated 
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Table 2: Noise components of CryoSTM-TIAs 


Pre-Amp type Single HEMT Amp Pre-Amp in Ref. [3] 
f (kHz) 10 100 10 100 
e? ((nV)*/Hz) 0.28 0.09 0.5 0.14 

Unit for the following terms is (fA)?/Hz 

in 0.14 4.9 0.8 20 
4kpT/ Rr 0.2 0.2 0.2 0.2 
B 0.0033 0.15 0.03 0.7 
(ae tei 0.28 0.09 0.5 0.14 
i2, as total 0.62 5.3 1.5 21 


in 


with the normal substitutional site configuration, called as shallow donor state; and a 
more localized acceptor level, known as DX~ center, arising from a lattice distortion. In 
Al,Gaji_zAs, as x > 0.22, DX” as a bound state is more stable than the shallow donor 
state. And, x = 0.37 for CNRS-HEMT in this work. There is a large repulsion barrier 
Ecap in the transition from the shallow donor to DX~, so DX? is frozen as T < 120 K 
27]. According to the work on the noise mechanism of the CNRS-HEMT [5], “bias- 
cooling method” can be used to control the density of the frozen DX- centers in the 
doping area in the HEMT. This unique method, by which the in-situ modification of the 
HEMT structure is realized, can be used to reduce the noise generated by the HEMT. By 
cooling the HEMT with the positive gate-source voltage, more electrons in the HEMT are 
frozen into the DX” centers. At low temperatures, the absolute value of the gate-source 
voltage required to maintain the HEMT at the ideal operating point is reduced, so that 
the gate leakage current is reduced, and thus the low frequency noises generated by the 
gate leakage current is reduced [5]. By the means, the inherent noise of the CryoSTM-TIA 
can be reduced. Moreover, if the HEMT used for the CryoSTM-TIA is not CNRS-HEMT, 
or its noise performances are not as good as those listed in Table [I] it can be remedied 
by “bias-cooling method”. This approach is very convenient for cases with only a single 
HEMT in Pre-Amp. 


4 Operating state adjustment and DC tunneling current 
measurements 


For the CryoSTM-TIA in Fig{I] along with the circuit parameters listed in Table [I] its 
operating state can be adjusted as follows. (1) Disconnect Pre-Amp from Post-Amp and 
ground the input of Pre-Amp (i.e. the Gate G of HEMT). Keep Rp to 0. Adjust Rs, Rn, 
and Rt, to achieve HEMT at the ideal operating point (Va, = 100 mV, Jg, = 1 mA) and 
the potentials at O; and Og equal. That is to say, Rg; and Ry, are adjusted to realize 
Ry = Rs + Vas/Ias = Rs +100 Q. (2) Cascade Pre-Amp and Post-Amp to form Inv-Amp, 
and input G is still grounded. Adjust Rp and Rs to realize the DC output voltage Vom 
of Inv-Amp to 0 and HEMT still at the ideal operating point. (3) Connect output O and 
input G of Inv-Amp with the feedback network and disconnect input G from ground, to 
form TIA. Input G is still not connected with the signal source circuit, so the potential 
of input G and output O of Inv-Amp is still 0, since the input resistance of the transistor 
HEMT Ra can be considered as infinity. In the feedback network, Rp + Rk + Rp = 1 
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GQ. (4) Connect the signal source circuit to TIA, to form the CryoSTM-TIA. 

When f — 0, aa > aao, and aao is the DC voltage gain of Inv-Amp. |aaol|an © 
(gmRf/g9aRun)ag © 94.5 dB, which is consistent with the simulation result |aao|an = 92 
dB shown in Fig[2| As the DC bias V; is applied by BMS, the DC resistance of TJ is 
R, and the potential at input G is Vg, and the output voltage of the CryoSTM-TIA is 
Vo. Obviously, aagVe = Vo, and (V; — Va)/R = (Ve — Vo)/Rr. The DC bias on TJ is 
V = V, — Vg, and the DC tunneling current I = (V; — Va )/R. J; as an approximate value 
of I is 

I, = —V,/ Re. (4.1) 


And, the relative error is obtained as 
Er = |I; — I|/|I| = 1/(1 — ago). 


As |aaolap © 92 dB, Er < 30 ppm, consistent with the simulation results [28]. Since 
aaoVa = Vo, the DC bias on TJ is 


V =Vi—Vo/aao. (4.2) 


By Eq. (4.1), Vo = -Rri ~ —Rel = —RpV/R. By Eq.(£2), Vi ~ V — ReV/(Raag). As 
the minimum of Ry is not less than 1 MQ, R > 1 MQ, so |Rp/(Raao)| < 1/40. Thus, 
V ~ V. With the measured V, I ~ J, and V ~ V; can be obtained. Therefore, the 
scanning tunneling current spectra J = I(V) (V €[Vi, Vu) can be obtained, where Vg 
and Vy are the lower limit and upper limit voltages for the measurements respectively. 

The drift of Vom can be considered as the amplification of the drift of the Inv-Amp 
input offset voltage Vos, where Vom = GagVos. For TIA, the output drift voltage is 
approxmate to Vog. The estimated drift of |Vos| is less than 17 wV/°C [29]. The power of 
the power supply of Pre-Amp is less than 40 mW. The temperature control system based 
on the TEC devices can be used to control the temperature fluctuations of the power 
supply and those of Rear-OPA within 0.01 °C, so the fluctuations of Vog within 170 nV 
is guaranteed, i.e. the TIA output fluctuations within 170 nV. 

In order to verify the correctness of the above design methods for TIA, we do the fol- 
lowing work shown in Ref.[3I]. The N-Channel JFET (SST4393-T1) is chosen to replace 
CNRS-HEMT, and TIA is designed and fabricated with the same circuit design method 
and the same circuit topology. SST4393-T1 can work at 77 K. We measured the electrical 
performances and noise performances of the JFET-based TIA at 77 K [BI]. The experi- 
mental results are basically consistent with the calculated results and simulated results, 
which verify the design methods for TIA in this work. 


5 Applications for the proposed CryoSTM-TIA in spectra 
measurements 

For most applications, the modulus of the transimpedance gain of the CryoSTM-TIA 

|A;(f)| should be measured firstly. Aj(f) is expressed as Eq.(2.7). Since |aa(f)|ap > 90 


dB in [1 kHz, 300 kHz] as shown in Fig] when Rj > 107° Rp, |Ai(f)| in [1 kHz, 300 kHz] 
can be approximately expressed as 


|Ai(f)| = Re/|1 — j2rfRp(Ca + Cr)/aa(f)|. (5.1) 
How to measure it has been described in Ref. [3]. 
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5.1 Measurements of scanning tunneling differential conductance spec- 
tra by the proposed CryoSTM-TIA 


The differential conductance of TJ Gy = 1/Ry is the function of the voltage V applied 
to TJ. As the frequency f of the modulated signal voltage Vi provided by BMS is low 
enough, such as f < 1 kHz, by Eq.(2.6), 

Rp 1 


Ay & . 
RV) 1- 


Rr i 
aa (f)R3(V) 


Therefore, Ry(V) ~ |-1/Ay + 1/aa(f)] Rr. When the measured |A,| < 1000, R3(V) ~ 
Rp /|Ay|, since |aa(f)|ag > 85 dB in f < 300 kHz as shown in Fig [2] With the measured 
|Ay| = |Vo|/|Vi|, the differential conductance spectra G3(V) = 1/R;(V) (V €[Vz, Va]) can 
be obtained. 

Increasing the frequency f of the modulated signal can speed up the scanning tunneling 
differential conductance spectra measurements. In [1 kHz, 300 kHz], as Ry > 1 MQ, by 
Eq. (2.6), 

(AIS po fir 

. IZa(f)| |1— g2afRr(Ca + Cr)/aa(f)| 
Therefore, 1/|Z3(f)| ~ |Av(f)|/|Ai(f)|, where | Ai(f)| is shown by Eq. (5-1). And, |Ay(f)| 
and |A;(f)| can be measured. 1/|Z3(f)| = 1/R3 + (2nfCy)* can be obtained. Selecting 
two different frequencies fı and fz in [1 kHz, 300 kHz], |Z3(f1)| and |Z3(f2)| are obtained. 
Ry and Cy can be solved out from the measured |Z 3(f1)| and |Z3(f2)|. 

Since the inherent noise of the CryoSTM-TIA is very small, the amplitude of the 
modulated signal voltage V, can be very small, so that the energy resolution for STS is 
much improved. 


5.2 Measurements of scanning tunneling shot noise spectra by the pro- 
posed CryoSTM-TIA 


The measurement method for the tunneling shot noise spectra by the proposed CryoSTM- 
TIA is basically the same as that introduced by Ref.[3], and it is only briefly described 
here. 

Before approaching the tip to the sample in the CryoSTM, Rj can be considered as 
infinity and Cy as 0. In this case, the output noise voltage PSD of the CryoSTM-TIA 
Ssu(f) can be measured, and the equivalent input noise current PSD i?(f) is 


P(f) = Sul PAA, (5.2) 
where |Aj;(f)| is obtained by Eq. (5.1). 


To measure the scanning tunneling shot noise spectra (STSNS) of a quantum system by 
the CryoSTM-TIA, the distance between the tip and sample is adjusted, and the interval 
(VL, Vu] is selected, so that the shot noise measurements at V € Dy = {V| V% < V < 
Vu,G3(V) < 1 uS} are needed to study the physical properties of the quantum system 
[3]. In Dy, the output noise voltage PSD Ssum(f, V) can be measured. The tunneling 
noise current PSD is denoted as S}(f,V). The equivalent input noise current PSD of the 
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CryoSTM-TIA is the function of Ry, so it is the function of V and denoted as i2 ( f,V). 
For Ssum( f, vV), 


Ssum(f,V) = [IEN EEN] AAP. (5.3) 
With Eq. (5.2) and (5.3), it is obtained that 


Ssum( f, V) _ Ssul f) 
|Ai(f)|? 


where 6(f,V) = 2,(f,V) — #2(f). From Eq. (3.12), 


Sil f, V) = o(f,V), 


R) =] + 


AkpT  e& = 
ae TA + (2r f)? (aa + 20CC A 


ateta 
Rr R ' 


And, i2 ( f, V) is composed by the four parts as shown in Eq.(3.12). According to Table |2| 
\5(f, V)| is smaller than 0.3 (fA)?/Hz in [10 kHz, 100 kHz) and 0.1 (fA)?/Hz in [100 kHz, 
200 kHz]. 6(f,V) can be neglected compared with S1(f, V), so long as the minimum of 


G3(V) (V € Dy) is not too small. Therefore, S1(f, V) can be obtained as 


Ssum aid — Ssu 
SIV) m 


by the measured Ssum(f, V), Ssu(f), and |Ai(f)|. And then, the SNTNS Sj,(J) = 2Fe|J| 
can be extracted from S{(f,V). Here, J = I(V) (V € Dy), and I is the DC tunneling 
current as the bias V is applied. 

Compared with the CryoSTM-TIA in Ref.[3], the CryoSTM-TIA proposed in this work 
has the same transimpedance gain, the bandwidth is above 200 kHz, and its inherent noise 
is much lower. The inherent noise of the CryoTM-TIA in this work is only 5.3 (fA)?/Hz 
at 100 kHz and its 6(f,V) < 0.1 (fA)?/Hz at 100 kHz, while that of the apparatus in 
Ref.[3] is 21 (fA)?/Hz at 100 kHz and 6(f,V) < 0.28 (fA)?/Hz at 100 kHz. Therefore, 
as investigating novel quantum states of various quantum systems, the measurements 
performed with the apparatus proposed in this work is more accurate. 

For example, as investigating the existence of MBS in a magnetic flux vortex of the 
iron superconductor in CryoSTM [4], the tunnel junction resistance must be large enough, 
since the occurrence of incoherent Andreev reflection can be ruled out by weak tunnel 
coupling conditions [32]. In Ref.[4], the tunneling current I is quite low (see Fig.3(a) in 
Ref.{4] and Fig.s2(a) and (b) in its supplemental file). As the tunnel junction bias V is 
0.3 mV, I is only tens of pA, and the corresponding shot noise may only be only a few 
(fA)?/Hz. For the shot noise measurements in the system in Ref.[4], it is obvious that the 
CryoSTM-TIA proposed in this work is much more effective and accurate than that in 
Ref. 3]. 


6 Conclusion 
In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning tunneling 


microscope (CryoSTM) is presented. TIA connected with the tip-sample component in 
CryoSTM is called as CryoSTM-TIA. The CryoSTM-TIA in this work has transimpedance 
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gain of 1 GO and bandwidth more than 200 kHz, and its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz. In this CryoSTM-TIA, a single CNRS-HEMT 
Pre-Amp is instead of the differential Pre-Amp consisting of a pair of CNRS-HEMTs in 
Ref.[B]. The difficulty of matching the identical HEMTs is avoided. And, due to the 
reduction of a noisy HEMT, the apparatus inherent noise is only 1/4 of that in Ref.[3}. 
Furthermore, the inherent noise of the single HEMT in the circuit can be in-situ reduced 
by “bias-cooling method”. With this apparatus, the fast high-energy-resolution scanning 
tunneling spectra measurements can be performed and the very low tunneling shot noise 
of quantum systems can be measured at the atomic scale. This apparatus can be applied 
to investigate novel physical properties of various quantum systems, such as detecting the 
existence of Majorana bound states in the topological quantum systems. 
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Supplemental file 1: Voltage gain peaks for cables 


Figure s1-1 The circuit scheme for measuring the voltage gain peaks for cables. 


A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance 
of 100 pF and a distributed inductance of 250 nH, which is shown as Fig. s1-1. In Fig. s1- 
1, R=1.5 kQ, the gain is Vo/Vin. The simulation results show its distributed LC structure 
results in several resonant gain peaks in [80 MHz, 1 GHz] shown as Fig. sl-2. As the 
input signal Vin is the thermal noise of R (1.5 kQ) at 300 K, we amplify Vo with 
Amplifier SA-220F5 (gain=200, fr=200 MHz), and then the gain peaks for the cable are 


measured with an oscilloscope. Fig. s1-3 show the measured gain peaks for cable. 
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Figure s1-2 The simulation results of the gain peaks for a 1 meter cable with a distributed capacitance 
of 100 pF and a distributed inductance of 250 nH. 
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Figure s1-3 The input signal Vin is the thermal noise of R (1.5 KQ) at 300 K. The measured results of 
the gain peaks for a 1 meter cable with a distributed capacitance of 100 pF and a distributed 
inductance of 250 nH. 


Supplemental file 2: Voltage gain of Inv-Amp 
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Figure s2-1 Inv-Amp circuit with all components for Eqs.(s2.1-s2.8), and their parameters are shown 
in Table 1 in Article. 


With the nodal analysis method, Inv-Amp voltage gain aa(f) is calculated by the 
following equations. All components for Eqs.(s2.1-s2.9) are shown in Fig.s2-1, and their 


parameters are shown in Table 1 in Article. 
(V, —V,)/R, =(V, —V,)/R, +(V, =E R +V,/Ry + joC,V, > (s2.1) 


JoCa (P-V) = 8a (P-V) + a (P-P) +0, -0 )/R +(V, -V.)/Ra +CV, 622 


( -P )/R +% -0 )/R, =V,/R, + j@C P> (82.3) 
(V,-V,)/R, =(V_-%)/R, JE s (s2.4) 
En (VV) + ga (P -8 )=P/R + JOCV, + jC, (V.-V,), (s2.5) 


V,=a,(V,-V), (s2.6) 


(V, ~ V.) /Ra =+joC,V, (s2.7) 


(V,-V.)/R, =+joC,V,, (s2.8) 


a, =V,/V, (s2.9) 
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Figure s2-2 The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|as, and the 
dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI simulation results, the green 
curves are the calculated results with the nodal analysis method (i.e. Eqs.(s2.1-2.9)), and the black 
curves are the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. 


Fig. s2-2 shows the modulus of the Inv-Amp voltage gain |aa(f)|az. The red solid 
curve is the TINA-TI simulation results, and the green solid curve is the calculated results 
with the nodal analysis method. They are identical with each other. The black solid curve 
is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In fe[3 kHz, 10 MHz], 


three curves are identical with one another. 


Fig. s2-2 also shows the argument of the Inv-Amp voltage gain Z(aa(f)). The red 
dashed curve is the TINA-TI simulation results, and the green dashed curve is the 
calculated results with the nodal analysis method. They are identical with each other. The 
black dashed curve is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In 


Jel kHz, 10 MHz], three curves are consistent with one another. 


Supplemental file 3: Transient response of the proposed CryoSTM-TIA 
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Figure s3-1 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with the 
different Ro for the time interval of 10 us. The dashed curve is the step input signal Vi/Vinst for the 
different Ro. The solid curves are the output response Vo/Vost. Here, Vost and Vinst are their values 
respectively at 300 ms after applying the step input signal Vi, which are shown in Fig.s3-2. Transient 
response time t<5 us. (b) is the zoom-in for the output response Vo/Vost in (a). 
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Figure s3-2 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 


simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with different 


Ro for the time interval of 300 ms. (a) The step input signal Vi/Vinst. (b) The output response Vo/Vost. 


Vo/Vos= 1 in t> 100 ms. (c) is the zoom-in figure of (b). 


A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. Fig.s3-1 shows TINA-TI simulation results for the transient response of the 
CryoSTM-TIA with the different Ro for the time interval of 10 us. The dashed curve is 
the step input signal Vi/Vins, and the solid curves are the output response Vo/Vost for the 
different Ro. Here, Vost and Vinst are their values respectively at 300 ms after applying the 
step input signal Vi. The error is definite as |Vo—Vost |/|Vost|. For the CryoSTM-TIA, the time 
taken from adding the input step signal voltage to the output response stably within a 
certain error is called transient response time f. As shown in Fig.s3-1 and Fig.s3- 
2(b)&(c), with Ro= 1 GQ and 30 MQ, for 0.1% error, t<5 us. With Ro=5 MQ, for 1% 
error, &<5 us. And, With Ro= 1 MQ, for 2.5% error, tr <5 us. 


Supplemental file 4: Noise of the proposed CryoSTM-TIA 


S4.1 Noises of STM-TIA 
For Inv-Amp, the equivalent input noise voltage of Inv-Amp is denoted as ea and its equivalent 


input noise current is ia, and their harmonic components of frequency f are Ea and Ja respectively. 


E 
(z can be obtained by the nodal analysis method. By Wiener-Khintchine theorem, 


A 


EE, Ex, 


ATA A` A 


2 a 
e ei 
A A^] can be obtained from ) [S4R1, S4R2]. The two matrix elements on 


ien tg 
the main diagonal are the equivalent input noise voltage PSD of Inv-Amp e} and its equivalent 


input noise current PSD i, . The two matrix elements on the sub-diagonal are its equivalent input 


noise voltage-current PSD e,i, and equivalent input noise current-voltage PSD i,e, . 


Inv-Amp is connected to the feedback resistor Rr to form TIA. The equivalent input noise 
voltage of TIA is denoted as er and its equivalent input noise current is ir, and their harmonic 
components of f are Er and Ir respectively. The temperature of the feedback resistance Rr is T, and 


its noise voltage is er and its harmonic component of fis Er. 


(b) 


Figure s4-1 (a) TIA circuit with the input short-circuit containing the equivalent input noise voltage of Inv-Amp 
ea and its equivalent input noise current ia, and the output noise voltage of eorv ; (b) Noiseless TIA circuit with the 
equivalent input noise voltage of TIA er as the input signal, and the output noise voltage of eotye; the equivalency 


of the above two circuits means €oty = €otve - 


For TIA, the circuit containing all noise sources with the input short-circuit is shown as 


Fig.s4-1(a), and the output noise is eoty. The noiseless circuit with the equivalent input noise voltage 
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of TIA er as the input signal is shown as Fig.s4-1(b), and the output noise is eotve. For calculating 
the equivalent input noise voltage of TIA, the equations are established on the equivalency of the 
above two circuits, i.€. Coty = @otve. Therefore, by nodal analysis method, 


ET=E,. 


Figure s4-2 (a) TIA circuit with the input open-circuit containing the equivalent input noise current of Inv-Amp 
Inv ea and its equivalent input noise current ia, and the output noise current of eori ; (b) Noiseless TIA circuit with 
the equivalent input noise current of TIA ir as the input signal, and the output noise voltage of eotie ; the 


equivalency of the above two circuits means eoti= €oric - 


For TIA, the circuit containing all noise sources with the input open-circuit is shown as 
Fig.s4-2(a), and the output noise is is eori. The noiseless circuit with the equivalent input noise 
current of TIA ir as the input signal is shown as Fig.s4-2(b), and the output noise is eotie. For 
calculating the equivalent input noise current of TIA, the equations are established on the 
equivalency of the above two circuits, i.€. eoTi = otic. Therefore, by the nodal analysis method, 


Tp=Iy+(Eat EpyVRe . 


(lun Jee) 


Hence, 


* * a 
From ee by Wiener-Khintchine theorem, m i can be obtained. The two 
LE, dhh ie i 


matrix elements on the main diagonal are the equivalent input noise voltage PSD of TIA e and its 


equivalent input noise current PSD ki The two matrix elements on the sub-diagonal are its 


equivalent input noise voltage-current PSD ei and equivalent input noise current-voltage PSD 


* 


iper: 
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oe i (s4.1) 


2 =i +4k,T/R,+e2/R2 , (s4.2) 


ei = =(re)* >ei +e fR, : (s4.3) 
Here, the noise voltage PSD of Rr is e = 4k,TR;. 
TIA is connected with the signal source circuit to form a STM-TIA. The equivalent input noise 


current PSD of a STM-TIA E can be obtained [S4R1-S4R3] by 


Pep +(1/R + +(2nf) Cher +(/R, + j2mfC, Jei +(I/R, — j2nfC, ire (s4.4) 
Here, Cu = CHC). Putting Eq.(s4.1), (s4.2), and (s4.3) into Eq.(s4.4), i? 


in =i, +4kgT/Ry +e% (YR +1/Re +(2nf) C3) 
(s4.5) 
+(I/R, + j2mfCy (enix +e} R; J+ (1/R, - j2mfCy) ine +e / Rs 


S4.2 Noises of the proposed CryoSTM-TIA 
S4.2.1 Equivalent input noise voltage and the equivalent input noise current of Inv-Amp 

For Inv-Amp, the equivalent input noise voltage and equivalent input noise current of HEMT 
are denoted as en and in respectively, and their harmonic components of f is are denoted as Ey and 
Iy respectively. gm is the transconductance of HEMT and ga is its channel conductance. The noise 
voltage of the resistors Ry is denoted as ern, and the harmonic components of fare denoted as Ern. 
The noise voltage of the resistors Ri, R2, and the feedback resistance Rr are denoted as ex1, er2 and ef 
respectively, and the harmonic components of f are denoted as Ei, Er2, and Er respectively. The 
noise voltage of the resistors Rc; and Rez are denoted as e1, e2 and ef respectively, and the harmonic 
components of f are denoted as E1, E2, and Ef respectively. Here, Ri= Ro= Rt. All noises of Ru, Ri, 
Ro, Rei, Re2, and Re above 1 kHz are thermal noise. Ra=Rt//(1/ga) and Rut = Rp//Rı are denoted. The 
equivalent input noise voltage and equivalent input noise current of Rear-OPA are denoted as ea and 
ia respectively, and the harmonic components of f are denoted as Ea and J, respectively. These noise 
sources are independent of each other. For Inv-Amp consistent of Pre-Amp and Post-Amp, its 
equivalent input noise voltage and equivalent input noise current as ea and ia, the corresponding 


harmonic components of fare Ea and Ja respectively. 
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Figure s4-3 (a) Inv-Amp equivalent differential circuit with the input short-circuit containing all noise sources, and the 
output noise current of eoay; (b) Noiseless Inv-Amp circuit with the equivalent input noise voltage of Inv-Amp ea as the 
input signal, and the output noise voltage of @oave; the equivalency of the above two circuits means eoay=oAve. In (a) 


and (b), the triangle OPA is the Rear-OPA. 


The equivalent differential circuit of Inv-Amp containing all noise sources with the input 
short-circuit is shown as Fig.s4-3(a), and the output noise is eoay. The noiseless circuit with the 
equivalent input noise voltage of Inv-Amp ea as the input signal is shown as Fig.s4-3(b), and the 
output noise is eoave. For calculating the equivalent input noise voltage of Inv-Amp, the equations 
are established on eoav = ove. Therefore, by the nodal analysis method, 

Ea ~ Eny + ErnZuv/(Avei Ru) EL (EmRL) + Za(E2/Ri + Ei/Re)/Avei 

+ Ea/Avpj + (E1-E2)/Ave1 ~ (1/8m + Zuu Avrj)ha , (s4.6) 
where Ayp1 = 2m/[1/Ritgaty(2nf/)Ci], Ave) =2m/[1/Ritgaty(2rf\Cir], Zar = Ruv/[1+ j2a)RarCir], and 
Cir =C; +. 

The equivalent differential circuit of Inv-Amp containing all noise sources with the input 
open-circuit is shown as Fig.s4-4(a), and the output noise is eoai. The noiseless circuit with the 


equivalent input noise current of Inv-Amp ia as the input signal is shown as Fig.s4-4(b), and the 
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output noise is eoaie. For calculating the equivalent input noise current of Inv-Amp, the equations are 
established on @oai= odie. Therefore, by the nodal analysis method, 
Ta~Tq t](20f)CarErnZur/(AveiRu) — J(27f/)CsaEui/(gmRi) + j2nf\CaiZar(Ev2/Ri+EdRs)/Aye 1 

+ j(2nf) (E1-E2)/Ayei— j(21f\CajE al Avr j (21) (Csa/gm + CajZu/A vej)la, (s4.7) 
where Csa = Cgs+Cga, Cat = CgstCga(1+Avp1) , and Caj = Cgst+Cga(1+Avp)). 


(a) 


Figure s4-4 (a) Inv-Amp equivalent differential circuit with the input open-circuit containing all noise sources, 
and the output noise voltage of eoai; (b) Noiseless Inv-Amp circuit with the equivalent input noise current of the 
Inv-Amp ia as the input signal, and the output noise voltage of eoaie; the equivalency of the above two circuits 


means €oai= odie. In (a) and (b), the triangle OPA is Rear-OPA. 


For the measured frequency of fm < 300 kHz, 2xfC;<0.3 mS with C=150 pF and 2zfC,<0.1 mS 
with C,=50 pF, but 1/Ra=ga+1/R1=2 mS and 1/Ryt=1/Ryt+1/Ri=6 mS, so Aypi=Pm/[1/Rit+gaty(2nf)Ci] 
and Aypj = gm/[1/Rit+gatj(2n/)Cir] are approximately simplified to —Ayp= gm/(1/Ri+ga) = 2mRa, Zui to 


Rut, and Caj to Ca. Thus, 
E, = 0 Top l Frc l Ry Een _ l Ey = l Rufo 
I) UJ* oy Wwe.) ARa ia )Cale.R (2a). j) ApR, (4.8) 
( l pes 4 l JAE -[ 1+ Ra /Ra j 
F(2mf Cy) ApRe UIO) C) Ar J (2af )(Ca + Ca Ra /Ra)) Sm 


By Wiener-Khintchine theorem, ignoring the small quantities, such as the thermal noise of R1, 
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Ro, and Rr, 


e ea H (Ri /Ra)em +(eter ter) |/45 +(14+ Ra /R) eles,» (s4.9) 

E R taal, Pa 

=i +(anp) co Hin Set + (rap RSE (rar)'[ C+ Fee, |, o 
Ri Ay Ap Ri m 


ae aN ee RÌ, ela etete R R P 
Tia en e (RERE eteo E] ee 
Here, E is the equivalent input noise voltage PSD of the CNRS-HEMT and 2 is their 


equivalent input noise current PSD. e? is the equivalent input noise voltage PSD of the Rear-OPA 


and T is its equivalent input noise current PSD. e? ; e2 , and ce are the thermal noise voltage 
PSD of Rei, Re2, and Rru respectively. 

e =e?=4k,T, R, =1.66 (nV)’/Hz , where Tr is 300 K. e=2.25(nV)/Hz and 
P=4 (pA)*/Hz in f> 10 kHz [S4R4]. In Eq.(s4.9)-(s4.11), (RR Jen is two orders of 


magnitude smaller than ¢?+e2+e2. (14R,, /R,) i?/e2} is one order of magnitude smaller than 
8 a e Te HL/ ta) la/Em g 


Z in Eq.(s4.9). (Ca+C, Ru /R, y E / gi is one order of magnitude smaller than 
Ci (e +e te, )/ 4 in Eq.(s4.10). (Ca+C, Ru /Ry)(1t+ Rap /R,)2/¢2, is one order of 
magnitude smaller than C, (e+e; +e] i A’, in Eq.(s4.11). Further ignoring the small quantities in 


Eq.(s4.9)-(s4.11), 


ee +(erte tez) /4& (s4.12) 
i pag +(2af) (oh (e +e? +e; e )/ A , (s4.13) 
ei = (ie) =- j(2af)C, (etere; A (s4.14) 


$4.2.2 Equivalent input noise voltage and the equivalent input noise current of of TIA 
For TIA, putting Eq. (s4.12) into Eq. (s4.1), putting Eqs. (s4.12) and (s4.13) into Eq. (s4.2), and 


putting Eqs. (s4.12) and (s4.14) into Eq. (s4.3), its equivalent input noise voltage PSD e2 and its 


equivalent input noise current PSD E , its equivalent input noise voltage-current PSD e,i% , and its 
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equivalent input noise current-voltage PSD i,e; are obtained as 


ete +(erte; +e | [4 (s4.15) 
iy = 4kpT /R, +e R? +(er ter +e; | NAR) +i + rY Cr lente: te;) /4. , (84.16) 
Cri =(i,€5) = 3 [Ry +(e e? i) (LR) -JOa)C (ete te) / 4 ' GFA 


54.2.3 The equivalent input noise current of the proposed CryoSTM-TIA 
Putting Eq.(s4.15), (s4.16), and (s4.17) into Eq.(s4.4), the equivalent input noise current PSD of 


the proposed CryoSTM-TIA A is obtained as 


2 =i +4k,7/R, +(1/R, +1/R;} lei (etete) [l +(2af) K ea +C (etete) Ja , 
(s4.18) 
where Cy = C1+Cy and C= Cat+Crt+Cj. Eq.(s4.18) is Eq.(3.11) in Article. And, 


i2 =i +4k T/R, +(1/R, +1/R f eè +i, , (s4.19) 


where i2, =(2nf) | Chek PEO (etete) i 4 | . Eq.(84.19) is Eq.(3.12) in Article. 


[S4R1] A. van der Ziel, Noise in Solid State Devices and Circuits, Wiley-Inter-Science, New York, 
(1986). 

[S4R2] Y.X. Liang, Low-noise large-bandwidth transimpedance amplifier for measuring scanning 
tunneling shot noise spectra in cryogenic STM and its applications, Ultramicroscopy, 234 (2022) 
13466. 

[S4R3] Z.H. Qian, Study on noise models algorithms and matrix descriptions for integrated circuits, J. 
Northeast Norm. Univ. 35 (2003) 41. 


[S4R4] Data sheet of THS4021 OPA, https://www.ti.com/lit/ds/symlink/ths4021.pdf.3. 
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Supplemental file 5: Estimating the DC tunneling current error for CryoSTM-TIA 
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Figure s5-1 For the proposed CryoSTM-TIA, Er = |L-I/|I| vs. R (Re[1 MQ, 1 GQ]) is simulated by TINA-TI with Vi 


=] mV (Black curve) and Vi=5 mV (Red curve). 


For the proposed CryoSTM-TIA, Fig. s5-1 shows the TINA-TI simulation results Er =|Js-|/|J| vs. R 
(Re[1 MQ, 1 GQ]) with Vi=1 mV (Black curve) and Vi=5 mV (Red curve). The simulation results 


show Er<30 ppm, consistent with the calculated results in Article. 
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Supplemental file 6: Estimation of Inv-Amp input offset voltage 


Figure s6-1 DC circuit of the inverting amplifier. 


Table s6-1 Parameters of THS4021 [S6R1] as Rear-OPA in Post-Amp 


Open loop DC voltage gain aao 97.5 dB 
Input offset voltage drift |vors/°C 15 V/C 
Input bias current Jop, Jon 3 uA 
Input bias current drift |ipp|/°C, |ibn|/°C Estimated as 100 nA/°C 
CMRR 95 dB 
dc = dao/CMRR 2.5 dB 


Pre-Amp in the left dashed box of Fig.s6-1 is made of CNRS-HEMT CH. Source of 


CH is grounded via a variable resistor Rs. Drain of CH is connected to Ri. A variable 


resistor Rp is adjustable from 0 to 10 Q. Ri=R2 = 


Ri =1 KQ. Its input (Gate of CH) is 


grounded. Its two outputs O; and O2 are disconnected from Post-Amp. PVR LM4050-10 


is used in the power supply to provide a constant voltage of 10 V (i.e. Vpp= 10 V here) 


with the typical temperature coefficient (T.C.) of 4 


t40) ppm/ °C [S6R2]. Keep Rp=0 and 


adjust Rr, Rs, and Ru to realize CH at the ideal operating point (Jas = 1 mA, Vas= 100 mV) 


and the DC voltage at O; equal to that at O2. The gate-source voltage Vgs =-—lasRs. For the 


ideal operating point, Vgs =—50 mV is assumed, so Rs =50 Q and Ru=150Q. 


Cascade Pre-Amp and Post-Amp to form Inv- 


Amp, and adjust Rr and Rp to realize 


CH still at the ideal operating point and the DC voltage at O to be 0. 


For Inv-Amp, assuming the output voltage drift of vo for 1 °C, the voltage drift at the 
output O; is vı, and it at Oj is v2. The voltage drift at the output of the power supply is vs, 
the voltage drift at PVR is vp, and the voltage drift at the source of CH is vs. By the 


following equations, the input offset voltage drift of Inv-Amp can be obtained. 


(vovi) RLH vv /RL= (Vp -)/Rrr , (s6.1) 
(vvi )/RL= ga(V1-Vs)-2mVs + ivp , (s6.2) 
(Vo-v2)/Ret (vi-V2)/RL= v2/Ru +ibn , (s6.3) 
Vs/ Rs = ga(V1-Vs) ~ZmVs , (s6.4) 
daol Vi tvo V2) + Ac(V1tVorr tv2)/2 = Vo , (s6.5) 
vos = Vo/AA0 . (s6.6) 
dao =~gmR/(gaRu) . (s6.7) 


where, vos is the input offset voltage drift of Inv-Amp for 1 °C. From the above equations, 
vos and vs can be solved out as vos ~ ibp(1+gmRs)/gm-inngdRu/gm-Vp(1+gmRs)/(2gmRr)- 
vof 1+2guRi+2mRs)/(2gmRi) and vs% vpRs/(2Rt)- ippRs. 

Therefore, |vp|<0.4 mV is the current drift of the power supply for 1 °C. High 
precision resistors with low T.C. are used in this work, so the drifts generated by the 
resistors can be neglected. Only considering the drifts generated by PVR and Rear-OPA, 
\vsi<10 uV, and |vos|<17 uV. The power of LM4050-10 is less than 40 mW in this work, 
and the power of THS4021 is less than 200 mW. It is easy to control the temperature 
fluctuations of LM4050-10 and that of THS4021 within 0.01 °C with the temperature 
Control System [S6R3], so the total input voltage drift of Inv-Amp can be controlled 
within 170 nV. 


[S6R1] Webpage of THS4021 OPA, https://www.ti.com/product/THS4021. 

[S6R2] Data sheet of LM4050-10, https://www.ti.com.cn/document-viewer/cn/lm4050- 
n/datasheet. 

[S6R3] A suitable TEC device, https://datasheets.maximintegrated.com/en/ds/MAX1978- 
MAX1979.pdf. 


2/2 


